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Abstract 



Using the Leiden/Dwingeloo Survey (Hartmann & Burton 1997 ^Uj) of the 
Galactic sky north of 5 = —30° as the principal component of a composite 
data cube, the structure of the warped HI layer of the outer Galaxy was dis- 
played by converting the data cube from heliocentric (I, b, -u)-coordinates to 
galactocentric (R, 9, z)-coordinates. We masked out known high-velocity- 
cloud complexes which might otherwise have contaminated the resulting 
description of the warped layer. We considered analogous displays under 
controlled circumstances by using as input medium a simulation of the same 
complete (I, b, v) data cube corresponding to a modeled HI Galaxy of known 
morphology and kinematics. By varying this artificial input we were able to 
check our method of converting the data cube and to put constraints on the 
global parameters that have been used to explain different asymmetries in 
the heliocentric data. 
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1 Introduction 

The early Leiden and Sydney surveys of HI in the Milky Way led to the dis- 
covery that the gas layer in the outer reaches of our Galaxy is systematically 
warped away from the flat plane 6 = 0°. Subsequently, it has become clear 
that many galactic gas disks are globally warped; indeed, it is difficult to 
argue that most galactic gas disks are not warped. The embedded perspec- 
tive which hinders some morphological investigations of our own Galaxy is 
in fact quite advantageous for studies of the warp. The first descriptions 
of the Milky Way warp as given by Burke (1957) 3], Kerr (1957) [12] : and 
Oort et al. (1958) [20] have not had to be altered in any fundamental way 
as investigations based on more modern data became available. The early 
work led to the recognition that the outer-Galaxy HI layer is warped above 
the equator b = 0° in the first and the second longitude quadrants, and 
below it in the third and fourth; the amplitude of the warp is larger in the 
northern data than in the southern; the total velocity extent in the south- 
ern outer-Galaxy layer is larger than in the northern; and the thickness of 
the gas layer in both hemispheres increases with increasing galactocentric 
distance. 

Although the essentials of the original findings remained unchanged, 
more recent studies by Burton & te Lintel Hekkert (1986) [Zj using prin- 
cipally the Leiden/Green Bank survey of Burton (1985) [S], and by Diplas 
& Savage (1991) (Hj using the Bell Labs horn antenna survey of Stark et al. 
(1992) [22], have extended the information and have yielded a more quan- 
titative description of the warp parameters. Also, several investigations 
have pointed out certain spectral asymmetries with respect to the gener- 
ally adopted approximation of kinematic and spatial circular symmetry to 
which a set of different explanations have been sought, not always in corre- 
spondence with each other (cf. Kerr (1962) [13], Blitz & Spergel, 1991 [2J; 
Kuijken, 1991 |E|; and Burton, 1991 0). 

The current project further extends the earlier work by using more sen- 
sitive data with more detailed spatial as well as kinematic coverage; it also 
attempts to account for possible contamination of the results by emission 
from high-velocity clouds. The data was transformed into galactocentric co- 
ordinates using the IAU standards of Ro=8.5 kpc and ©o=220 km/s, and a 
wide range of morphological parameters was further quantified. Comparing 
observationally established spectral asymmetries with respect to kinematic 
and spatial circular symmetry with similar asymmetries in synthetic spectra 
for a range of Galactic morphological and kinematic situations, we attempt 
to identify the response of various aspects of the spectra to small changes 
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Figure 1: Representative cuts through the all-sky composite datacube. The 
upper panel shows an (l,v) plot at b = —2°; the lower panel, at b = +2°. 
The warp at negative latitudes in the third and fourth quadrants causes the 
emission to extend to larger velocities when probed below the plane 6 = 0°, 
and vice versa for the warp at positive latitudes in the first and second 
quadrants. The peak intensity corresponds to a brightness temperature 
around 135 K. The lightest greys in these displays correspond to Xft=4 K, 
much greater than the sensitivity of the LDS, cr=0.07 K. All (I, v) plots in 
this paper use the same grey-scale, except for those latitude-averaged. 

in these models and to put constraints on the parameters used in them. 

2 Data used 

The project utilizes the Leiden/Dwingeloo HI survey of Hartmann & Bur- 
ton (1997) which has a greater latitude coverage than that of the Lei- 
den/Green Bank survey (which was restricted to \b\ ~ 20°), and a greater 
sensitivity (a = 0.07K); the new data have a more detailed spatial cover- 
age (AZ = 30', Ab = 30') as well as a more detailed kinematic resolution 
(Av = 1 km/s) than the Bell Labs data and were corrected for contamina- 
tion by stray-radiation. The Leiden/Dwingeloo data covers the sky north 
of 5 = —30°. HI data of the full southern sky are not yet available, although 
the latitude range |6| < 10° at 5 < —30° is covered by the Maryland/Parkes 
survey of Kerr et al. (1986) |14| : in the region below 5 = —30° where 
the survey of Liszt & Burton (1983) was available, it was given prece- 
dence over the Parkes data because of its higher resolution. The data from 
the two additional sources were interpolated onto the (l,b,v) grid of the 
Leiden/Dwingeloo survey and scaled to brightness temperature using the 
parameters of the S7 standard field given by Williams (1973) j2&| . In Figure 
1 the resulting all-sky data cube was displayed for two cuts through galactic 
longitude, I, and velocity, v, at constant galactic latitudes of b = —2° and 
b = 2°. 

3 Masking high-velocity-clouds 

The distortions by high-velocity-cloud material were addressed by masking 
out the principal HVC complexes. We note that in some directions at inter- 
mediate latitudes the total HI emissivity contributed by anomalous-velocity 
gas associated with the high-velocity and intermediate-velocity clouds can 
be greater than the conventional- velocity emissivity. Unless this material 
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Figure 2: Rotational situation used in the circularly symmetric velocity-to- 
distance conversion of the data. 

is masked out, the velocity-to-distance analysis in section four would have 
attributed it to distances calculated according to the normal rotation veloc- 
ity (unless the HVC velocity were so extreme as to be formally forbidden 
in terms of the input rotation). To test what the effect of this disturbance 
is on the displays of the data after conversion to galactocentric (i?, 9, z) 
coordinates, we masked out the major HVC complexes A, C, H, and AC, 
as well as the complexes WA, WB, WC, and WD (see e.g. Wakker & van 
Woerden 1997) |24j . Complex G was partially left behind in order not to 
loose any well-behaved data since it has low deviation velocities. Then the 
masked data set was subtracted from the observed data and the result was 
displayed in the (R,0,z) coordinates. The differences between the results 
obtained without masking out the high-velocity clouds and those obtained 
after masking were minimal. The velocities of the high-velocity cloud com- 
plexes are largely forbidden by any reasonable rotation curve. Therefore 
we conclude that the properties of the Figure 3-6 type of displays are not 
significantly affected by HVC contamination. 

4 Coordinate conversion 

The brightness temperatures in the masked data cube of heliocentric co- 
ordinates (l,b,v) were then transformed to HI volume densities in galac- 
tocentric cylindrical coordinates (R, 9, z) using a kinematic model of the 
Galaxy that is circularly symmetric and has a flat rotation curve with 
0(i?>R o =8.5 kpc)=220 km/s. Figure 2 shows a face-on view from the 
northern side of the Galaxy of the rotational situation. An observer at the 
sun's position, embedded in the equatorial plane of the rotating Galaxy, 
measures a radial velocity along its line of sight towards an object at po- 
sition (l,b,R) of: v(l,b,R)=R sml[@(R)/R - e(R )/R ]cos b. After the 
velocity-to-distance transformation was carried out, volume densities were 
derived in galactocentric cylinders, assuming optical thinness. Since at any 
specified velocity range the optical depth, t(v), and the column density, 
Nhi(v) (in number of atoms at a particular velocity in a cylinder of cross- 
sectional area 1 cm 2 ), are related by N HI (Av) = 1.823 x 10 18 / T k T(v)dv 
[cm -2 ], with Tfc the spin temperature of the gas. For the case of optical 
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thinness, T b (v) = T k (l - exp- T ^) ~ T k r(v), so we can derive the volume 
density through the integrated intensity observed: tihj(Av) oc fT{,{y)dv/Ar 
[cm~ 3 ]. If, on the other hand, the HI in some direction contributing at a 
particular velocity is opaque (t(v) S> 1), the observed brightness temper- 
ature does not reveal the number of emitting atoms, but, instead, only 
a lower limit. It does indicate, however, the gas temperature, because 
Tb(v) = 7fc(l — exp" 1- ^) ~ Tfc. In this case, the profile integral yields an 
interesting measure: / Tf,(v)dv J Ar ~ T k \Av / Ar\. The expression \dv/Ar\ is 
called the velocity-crowding parameter. It is a measure of how fast the line 
of sight distance varies for a given radial velocity increment, according to 
the pertaining kinematics. Emission from regions of high velocity crowding, 
that is for a low absolute value of dv/Ar, is compressed into a short veloc- 
ity interval making the spectra tend to saturate. This effect is particularly 
notable in both cardinal directions (I = 0° and I = 180°) and also, though 
somewhat less, in the directions I = 90° respectively / = 270°. 

The volume densities resulting from the conversion are shown in Figures 
3 through 6 in cuts through azimuth, 9, (defined such that 9=1 for a point at 
infinity) and vertical height, z, at constant galactocentric radius, R, at az- 
imuth intervals AO = 1°, in vertical thickness steps Az=25 pc, in cylindrical 
walls of thickness Ai?=250 pc. 

We restricted the use of data to those latitudes that are essential in 
order to display the warp to any sensible extent in radius. Thus we used 
data with b > —20° and b < 30°. This latitude cutoff can be seen clearly 
in the displays as arcs through z and 9 at smaller radii. One immediately 
also sees the limited extent of the Maryland/Par kes survey of Kerr et al. 
(1986) [Hj (\b\ < 10°) as well as that of the Liszt & Burton (1983) [TH] survey 
(|6| < 10°, with some indentations at the most negative latitudes). Figures 
3 through 6 show that these latitude limitations put no major restrictions 
on the validity of the derived warp parameters, as the data in the northern 
hemisphere up to b = 30° seems quite sufficient to display the warp and, 
due to the warp's lesser amplitude in the southern hemisphere, so does the 
data down to b = —10°. 

At 9 = 0° the emission shows a dip. This is not an observed feature, 
but a consequence of the conversion's analytics that reach a mathematical 
singularity at 9 = 0°. Also, although the displays still show non-zero HI 
densities at radii larger than the last displayed at R = 26 kpc, they are not 
shown here because of large systematic errors. At those radii the densities 
have lessened to such extent that it is difficult to distinguish between a 
density decrement with increasing radius as a galactic morphological artefact 
and a density decrement along velocity dispersion wings which is converted 
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Figure 3: Mean HI volume densities in cylindrical cuts through the compos- 
ite and HVC-masked data cube at distances of R = 12 kpc (upper) and 14 
kpc (lower). The cylindrical-coordinate (R, 9, z) data set was sampled at 1° 
intervals in 6 and 25-pc intervals in z; the velocity-to-distance transforma- 
tion displayed here was circularly symmetric. The darkest grey-scale was 
scaled to the maximum densities at each radius, to be able to show details 
even at large R, and can be compared with uhi values in Figures 9 and 10. 

Figure 4: Same cuts as in Figure 3, at galactocentric distances of R = 16 
kpc and 18 kpc. 

as a progressively smaller amount of gas originating from larger radii (Figure 
12). It is further evident in the displays that the heliocentric sampling of 
space reaches finer gridpoints when getting closer to the sun by noticing that 
a higher resolution is attained closer towards the anticenter. 

In order to follow the transformation process under controlled circum- 
stances, we have simulated the complete observed (l,b,v) data cube on the 
same grid by solving the radiative transfer pertaining for an ad hoc Galaxy 
characterized by a plane-parallel, but warped and flaring, gas layer. This 
both provides an objective control on the conversion method used and gives 
at the same time useful insight in the disturbing effects on the images in 
(R, 0, z) coordinates of velocity crowding, changing degrees of optical thin- 
ness, and a finite velocity dispersion of the gas. More on this will follow in 
sections 5.2.1 and 5.2.2. 

4.1 Derived galactic morphology 

After conversion, the resulting set of galactocentric representations of the 
Galaxy, assuming circular symmetry, can be used to obtain a first-order 
impression of its morphology. For comparison with earlier data we computed 
several parameters concerning this morphology. At the same time these 
parameters can be used as the most realistic input values currently available 
for the simulations of section five. It must be noted, however, that an 
unavoidable error is being made here inherent to the method, when we make 
statements on possible galactic morphologies and kinematics as responsible 

Figure 5: Same cuts as in Figure 3, at galactocentric distances of R = 20 
kpc and 22 kpc. 
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Figure 6: Same cuts as in Figure 3, at galactocentric distances of R = 24 
kpc and 26 kpc. 

for observed asymmetries on the basis of simulated models that took as input 
the parameters that were computed assuming circular symmetry. Since all 
parameters were azimuthally averaged to fit circular data and, moreover, 
since even in a worst case scenario the errors will typically not exceed a few 
to at most ten percent, we assume that this has no important effect on the 
conclusions that follow. 

The parameters in the Figures 8 through 12, and 14 and 15 were com- 
puted by fitting a Gaussian to the z-dependent HI density distribution for 
each fixed (R, 9) combination, using the GAUSSFIT routine in IDL. We 
measured for each (R, 9) combination the height of the Gaussian, its width 
and the location of its center. From the z location of its center we deter- 
mined the parameter plotted in Figures 8 and 9, from its height, the one 
in Figures 10, 11, and 12, and from its width, the one in Figure 14 and in 
Figure 15. It should be noted that the Gaussian width corresponds to the 
dispersion a and not to a full width at half maximum (FWHM=2.36*<r). 
The parameter in Figure 13 was computed by integrating all z values for 
each point in (R, 9). 

It would be valuable, of course, to be able to determine the total extent 
of the Milky Way's HI disk in view of its theoretical consequences for, e.g., 
the shape of the Galactic dark halo and the Local Group's dynamics, to 
name but a few. However, this is not easily established. Not only is there 
the difficulty to distinguish in the observed spectra between the outer edge 
of the Galaxy and its less outward velocity dispersed gas (of which the exact 
dispersion rate at these outer parts is difficult to determine) as is described 
above, but also is it likely that errors in this determination are superseded by 
errors in the assumption of the Galactic kinematics as used in the conversion 
method. It is clearly of influence what type of rotation curve one uses in 
the velocity-to-distance transformation; currently, there is not one uniquely 
determined rotation curve found yet. Also, there remains some uncertainty 
in the correct value found for Ro Furthermore, as will be commented upon 
in section 6.1 and 7.7, some deviations from kinematic and spatial symmetry 
have important consequences for the perceived Galactic morphology after 
conversion. In particular, note that it appears that the HI density at large 
radii is large in the fourth quadrant with respect to the first and that of the 
second quadrant is large compared to the third. As will be shown (in section 
7.7) this effect is also reproduced with slightly elliptical HI streamlines in 
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Figure 7: Line of sight velocities for emission from increasing radius R in the 
direction I = 90° ,6 = 0°. The solid line is according to the kinematics used in 
this investigation (Q(i?>R =8.5 kpc)=220 km/s), the dotted line to that in 
Henderson et al. (1982) [TT] and Diplas & Savage (1991) (G(i?>R o =10.0 
kpc)=250 km/s), and the dashed line corresponds to the kinematics specified 
by R =8.5 kpc, G o =220 km/s, and Q(R) = G ( J R/R ) 0382 as is used by 
Wouterloot et al. (1990) |H] 

Figure 8: Vertical distance from the plane b = 0° of the gauss-fitted peak HI 
density at galactocentric distances of respectively R = 12 kpc, 16 kpc, and 
20 kpc. The distribution was smoothed over 10° in 9. The central bump 
near 9 = 0° is not physically significant, due to the smearing of the data in 
that direction (see the comments on this in the caption of Figure 14). 

Galactic simulated spectra. That the data, and the parameters deduced 
from it, are nevertheless displayed in (R,6,z) coordinates, using cylindrical 
symmetry and the rotation curve as above, is for comparison with earlier 
data (to which this has the advantages described in the introduction) and by 
lack of any better established alternative yet. The same arguments lead to 
the explanation of the absence of error bars on the parameters displayed here, 
as it is estimated that the systematic errors following from the conversion 
method are greater than the errors introduced by fitting Gaussians or than 
errors in the data itself. 

It is instructive to compare these new data with the findings of past 
investigations. For instance, although the shapes of the surface density (a) 
curves of Henderson et al. (1982) in their Figure 4 and that of our 
Figure 13 are quite similar it can be seen that the peak of their northern 
surface density profile occurs at a radius of ~14 kpc, while ours occurs at 
~12.5 kpc. This difference is caused by the choice of the assumed Galactic 
kinematics in the velocity-to-distance transformation as is illustrated in 
Figure 7. This scale effect, which is more pronounced for radii nearby to the 
solar radius, is squared when it enters their analysis of the surface density, 
using G(i?>R o =10.0 kpc)=250 km/s, and leads to an underestimate relative 

Figure 9: Vertical distance from the plane b = 0° of the gauss-fitted peak 
HI density at galactocentric distances of respectively R = 24 kpc, 26 kpc, 
and 28 kpc. The distribution was smoothed over 10° in 9. 
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Figure 10: Azimuthal distribution of HI midplane densities, measured as the 
gauss-fitted peak of the vertical HI density distribution, at the respective 
radii of R = 10 kpc, 12 kpc, and 14 kpc. The distribution was smoothed 
over 10° in 9. The plots show an apparent underabundance of HI gas in 
those directions (9 ~ 0°,±180°) where due to velocity crowding the total 
amount of emitting gas is expected to be higher than what is observed by 
received emission. Note, specifically, the low values around / = 0°. 

Figure 11: Same azimuthal distribution of HI midplane densities as in Figure 
10, at the radii of R = 16 kpc, 18 kpc, and 20 kpc. 

Figure 12: Azimuthally averaged HI midplane densities with increasing 
radius for the range 20° < 9 < 160° (dotted line) and for the range 
200° < 9 < 340° (dashed line). It is difficult, at very large radii, to dis- 
tinguish between a density decrement with increasing radius as a galactic 
morphological artefact and a density decrement along velocity dispersion 
wings which is by the conversion method naturally taken as a progressively 
smaller amount of gas originating from larger radii. 

Figure 13: Azimuthally averaged HI surface densities (a) in solar masses 
per square parsec (upper panel), for the range 20° < 9 < 160° (dotted line) 
and for the range 200° < 9 < 340° (dashed line). The lower panel gives the 
logarithmic version from which we deduce that the mean radial scalelength 
for the HI surface density in the 6*-averaged range 20° < 9 < 160° at R > 9 
kpc is 4.3 kpc and for the range 200° < 9 < 340° 4.9 kpc. However if we 
only measure at 15 < R < 25 kpc, then we arrive at 3.7 kpc and 4.2 kpc, 
respectively. In view of the low densities at R > 25 kpc, it is difficult to 
determine the significance of the values of a at those radii. 

Figure 14: Gaussian dispersion (h) of the vertical HI density distribution 
at distances of R = 12 kpc, 16 kpc, and 20 kpc. The distributions were 
smoothed over 10° in 9. Note the increase of average scale height with 
increasing radius. The large values of h near 9 = 0° are spurious, being due 
to the elongated patch of emission in the Figures 3-6. This is caused by 
local gas that reaches velocities, due to the finite velocity dispersion, that 
correspond kinematically with radii in the far outer Galaxy. Note, that it is 
slightly off-center, as is the elongated patch. 



4 COORDINATE CONVERSION 



12 



Figure 15: Flaring of the azimuthally averaged gaussian width (h) of the 
vertical HI density distribution for increasing galactocentric radius. The 
dotted line is the ^-averaged value for the range 20° < 9 < 160°, the dashed 
line is the ^-averaged value for the range 200° < 9 < 340°. Finally, the 
solid line is their average. The 40°-bands in both cardinal directions were 
omitted, because of their unrealistic high values of the scale height as ex- 
plained in the caption of Figure 14. As mentioned in Figure 13, also here it 
is difficult to determine the significance of the values of h at R > 25 kpc, in 
view of the very low densities at those radii in Figure 12. 

to the analysis we made, using 0(i?>R o =8.5 kpc)=220 km/s. Thus the 
difference in the magnitude of the density, with the current findings coming 
to a more than 2 solar masses per square parsec higher value than those of 
Henderson et al. for radii up to about 17 kpc, is at least in part accounted 
for. This is clearly also the case in the logarithmic plot of a in Figure 8 of 
Diplas & Savage (1991) who used the same kinematics as Henderson et 
al. Also, they show in their table 1 that derived scale lenghts are expected 
to be greater when using the kinematics we used than using theirs. They 
point out that for the Henderson et al. data they estimate r s =2.08 kpc 
(for 9( J R>R o =10.0 kpc)=250 km/s), while they arrive at r s =2.49, 2.60, 
and 3.98 kpc for values of 9 of 50°, 90° respectively 130° corresponding to 
r s =3.14, 2.74, and 5.55 kpc for 9( J R>R =8.5 kpc)=220 km/s. These three 
specific values can not be compared directly to the ^-averaged values we 
found as explained in the caption of Figure 13, though it is clear that ours 
are significantly higher than the one found from the Henderson et al. data. 
Since the shape, and particularly the height, of the Galactic surface density 
curve is very sensitive to the choice of range for averaging in 9, as can clearly 
been seen in Figures 10 and 11. In the case of the work done by Henderson et 
al. this range was 30° < 9 < 150° (northern) respectively 210° < 9 < 330° 
(southern) where we used 20° < 9 < 160° respectively 200° < 9 < 340°, 
some of the difference in surface densities along Galactic radius may also be 
accounted for by this effect, although the lower densities toward the cardinal 
directions might in fact argue for a bias in the Henderson et al. data towards 
higher values with respect to ours, instead of lower. That this effect is 
significant and can be attributed to velocity crowding alone is shown in the 
lower panel of Figure 22 in section 5.2. This is also the direct reason for us 
to leave out a band of 40° in 9 in both center and anticenter direction in 
our analysis. The data as analyzed by Wouterloot et al. (1990) |26j is more 
directly comparable to the current results, as can be found in their Figure 5b, 
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due to the small difference in kinematic distance determination between their 
method, using R Q =8.5 kpc, 6 o =220 km/s, and &(R) = o (i?/R o ) a0382 , and 
ours, as can again be seen in Figure 7. They find for the ranges 120° < 9 < 
170° (northern) respectively 190° < 9 < 240° (southern) equivalent, though 
slightly lower, magnitudes than ours, where due to a bigger overlap with the 
region that is most affected by underestimation of density through velocity 
crowding one would expect these magnitudes to be lower limits, rather than 
upper. This is almost certainly the case for the Galactic surface density curve 
shown in their Figure 5e, where they averaged over the complete second 
and third quadrants and thus included the minimum around / = 180°. As 
therefore expected, it displays lower a, especially at radii nearer to Ro=8.5 
kpc, which can be understood when taking into account the progressively 
larger relative minimum in the Figures 10 and 11 for values of R nearer to 
the solar radius. Their value of radial scalelength, r s = 4.0 kpc, is, again, 
equivalent to ours, though slightly lower, in spite of the suppressed density 
magnitudes at radii close to Ro, which would make it larger. This could 
be attributable to the fact that at larger radii the difference between the 
kinematics becomes more pronounced (Figure 7), causing their a to decline 
faster. 

The flaring of the thickness of the HI layer is demonstrated in Figure 
15. Also here it can be seen, from Figure 14, that the exact shape and, in 
particular, the height of the flaring curve is sensitive to the applied range 
of ^-averaging. It is remarkable how the data from Wouterloot et al., in 
their Figure 9, leads to a much smaller thickness near R = 10 kpc of a 
half width to half maximum of about 200 pc, which roughly equals 170 pc 
(FWHM=2.36*cr) in Gaussian width (a), compared to our ~450 pc in a. 
However, Diplas & Savage find, using the data of Henderson et al., for a 
(9-averaging over 90° < 9 < 130° about cr=500 pc and for 210° < 9 < 250° 
about cr=400 pc, both at R = 12 kpc, in good agreement with the data 
here. Their plots of the flaring of the thickness for different particular 9 
values in this figure, with their wide range in thickness at any particular 
radius, illustrate the impact of the range of 9 included on the flaring curve, 
as was mentioned above. We choose not to include a band of 40° in 9 in 
both center and anticenter direction in our analysis, since we can see that 
in these regions the values of h are unrealistically high due local velocity 
dispersed gas, explained in the caption of Figure 14, as is demonstrated in 
Figure 22, upper panel, in section 5.2. For all the plots showing the flaring 
of the thickness, the flaring continues to rise towards higher radii. Only 
our new data shows a slowing of rise in flaring at very large radii (R > 23 
kpc), although one should be careful about the interpretation of the analysis 
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Figure 16: Solid lines trace the sine-fitted warp heights (see text) for the first 
and second quadrant together (upper) respectively third and fourth (lower). 
The dashed lines are the mathematical fits that are used in the model. 

at these large distances from the Galactic center in view of what was said 
before about systematic errors. 

When looking at the vertical distances from the plane b = 0° of the 
gauss-fitted peak HI density at very large galactocentric distances in Figure 
9 the attention is drawn to a feature, which is referred to in the literature 
as scalloping (e.g. Kulkarni et al. (1982) JZ1); that is most pronounced in 
the 3 rd and 4 th quadrants. As mentioned before, carefulness is required for 
the interpretation of the reality-content of this feature at such large radii. 
Closer inspection reveals that about 10 periods can be distinguished for a 
full 360° in azimuth, the same number that was found by Kulkarni et al. 
The scalloping is further remarkable, because it sustains the same regularity 
through several kpc in galactocentric radius. Currently, it is uncertain, what 
its cause is. In any case this effect is additional to the global parameters we 
consider here and shall not be pursued by us any further. 

4.2 Computed input parameters 

We have constructed a complete set of (l,b,v) spectra based upon an HI 
brightness temperature distribution according to known input galactic mor- 
phology and kinematics. More on this will follow in section 5. Both for the 
purpose of having a realistic controlled input medium for the heliocentric- 
to-galactocentric coordinates conversion and for the purpose of showing ob- 
servational consequences of various different galactic models that have been 
put forward as a possible explanation of observed spectral asymmetries, it is 
valuable to have close-to-reality input parameters. Some of these are taken 
from the literature, others have been computed by using this project's data 
and the parameters calculated in the subsection above. We will mention 
those of the last type here and explain the method by which they were 
calculated. 

The radial scale length was taken as an average of the four values that 
were found from the lower panel in Figure 13. It was derived to be 4.3 
kpc. The vertical scale height in the solar neighborhood was deduced from 
Figure 15. This is a Gaussian scale factor and it was computed to be 0.40 
kpc. The rate of flaring was also taken, azimuthially averaged, from the 
same plot and amounted to 25 pc in Gaussian width of vertical thickness 
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for every kpc radially outward beyond R =9 kpc, the radial edge of the 
assumed flat and of constant density disk. This central disk HI density was 
measured at a radius of R = 9 kpc and averaged over all 9 values except for 
20° regions centered on both the anticenter and center directions in order 
to avoid an underestimation due to velocity crowding. We found 0.26 cm~ 3 , 
which is a lower limit. At a certain point the simulation makes use of a 
sinusoidal warp that is linearly rising beyond R = 9 kpc. That this is a 
fair first approximation can be seen in Figure 8. To determine the rate of 
rising of the warp we calculated the magnitude of a sine curve that maps 
out the same surface in a plot either above or below the plane 6 = 0°. These 
thus determined sine-fitted warp heights for increasing radius can be found 
as the solid lines of Figure 16. Their increasing rate was averaged to give 
110 pc per kpc radially outward from R = 9 kpc. A better approximation 
involves an asymmetric warp that distinguishes a part at latitudes north 
of the plane 6 = 0° and one south of it. This approximation increases not 
linearly but according to the mathematical fits shown with dashed lines in 
Figure 16. 

In the original (I, b, v) data the warp maximum and warp minimum were 
defined as the longitudes at which the highest respectively lowest latitudes 
were attained. This was measured for a velocity-integrated (l,b) plot at 
—220 < v < —100 km/s for the warp at positive latitudes in the first and 
second quadrant to include only gas from the outer Galaxy, as well as for 
100 < v < 220 km/s for the warp at negative latitudes in the remaining two 
quadrants. This resulted in / = 103° for the warp maximum and I = 262° 
for the warp minimum. Because of the smaller than 180° spacing between 
the maximum and the minimum and the known folding-back property of 
the warp at the southern hemisphere, which is prominent in the velocity 
range that was sampled here, the measurement was repeated for the range 
60 < v < 100 km/s, which is representative for distances more towards those 
where the southern warp reaches its minimum, shown in Figure 16. This 
gave for the warp at b < 0° a quite shallow minimum that extended from 
/ ~ 240° to I ~ 290°. A warp maximum for / = 100° implies for the models a 
galactocentric warp maximum at 9 ~ 120° spaced from the warp minimum 
by 180°, assuming that this maximum is observed at R ~ 20 kpc or more. 
However, judging by the line of nodes in Figures 8 and 9 one would come 
to a warp maximum at 9 ~ 95°. As, from the same figures, it is clear that 
both orientations represent only in part the real situation, since the shape 
of the warp is not perfectly sinusoidal, both will be used in the simulations 
of section five and seven. 
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Figure 17: General elliptical rotational velocities used in the simulation 
models, Q x has a negative value, Q y a positive one. Note that in the case 
e = and a = 0°, it will follow that 4> = 0. 

5 Modeling under controlled circumstances 

We have constructed a complete set of (l,b,v) spectra based upon an HI 
brightness temperature distribution according to known input Galactic mor- 
phology and kinematics. This simulation provides a strong check on the 
analysis procedures, when taken with identical kinematics as those specified 
in the method used for converting the data to (R, 0, z) coordinates. Chang- 
ing the input models, it is furthermore possible to monitor the influences 
this has on both the appearance of the spectra, and on the appearance of 
plots in (R, 8, z), that were made, while still using the conversion we use in 
the real data's case. To be able to include demonstrations of the effects of 
any possible kinematic ellipticity or other forms of asymmetry, we have used 
as an input model the general case of a Galaxy with kinematics described 
by motions in closed elliptical streamlines. That it is reasonable to do so, 
even though in a non-spherical potential there are no dynamical stable and 
closed orbits, can be thought of when one realizes that a collection of HI 
clouds on independent and randomly chosen non-stable orbits are likely to 
collide and thus to cause shocks, for which there is no observational evi- 
dence. The kinematics were therefore chosen without a formal dynamical 
foundation, with an ellipticity, e, defined for the rotational velocity, 0, of 
any object at a certain position angle ip (which is a + 180° in the case of the 
sun; see Figure 17). This is decomposed into a velocity along the x-axis 
(defined by the ellipse major axis), Q x , and a velocity along the y-axis, Q y , 
where Q x = — 0o(l + e/2) sin tp and 0^ = ©o(l — e/2) cos if>. 

5.1 Circular symmetry 
5.1.1 Description of the model 

The simplest case we consider is that of a plane-parallel circularly symmetric 
Galaxy, with differential rotation according to @(R) = 0o[l.OO74(i?/R o ) ao382 + 
0.00698]-2.0 (km/s), for R <R , and @(R) = O , for R >R . The gas layer 
shows flaring, but is not warped. The midplane is specified by the plane 

6 = 0° and has a volume density of 0.30 cm -3 , a little higher than what 
was found in section 4.2 to compensate for the underestimation of this pa- 
rameter due to the fact that not everywhere in the Galaxy along the line 
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Figure 18: Simulated analogue at b = 0° of Figure 23, upper panel, for 
a model with complete circular symmetry as described in the text. As 
in Figure 1, the peak intensity of the display corresponds to a brightness 
temperature of 135 K; the lightest greys, to Tfe=4 K. 

of sight the opacity is truly far smaller than 1. That this expectation is 
justified can be seen in the middle panel plot of Figure 21, which bears a 
fair resemblance to that of Figure 12. The volume density is constant at the 
midplane up to R = 9 kpc, but at larger R it declines with a scale length of 
r s =4.3 kpc. Moreover, it also declines with vertical \z\ height, with a scale 
height of /i s =0.40 kpc. Due to the flaring of this thickness, this value of 
h s increases with 0.025 kpc for every kpc in radius beyond R = 9 kpc. All 
these parameters were derived using the real data after a circularly symmet- 
ric conversion as described in section 4.2. The model Galaxy is truncated 
at R = 30 kpc. The gas was defined to have a velocity dispersion of a = 7.5 
km/s and a spin temperature of Tj, = 135 K. 

5.1.2 Simulation results 

The complete observed (l,b,v) data cube is then simulated on the same 
grid by solving the radiative transfer pertaining for the ad hoc Galaxy, as 
described above. Figure 18 shows the resulting set of spectra at b = 0° in 
a display through I and v. Note that, apart from the scale factors men- 
tioned above, the input model did not have any structure and that thus all 
artefacts in this (I, v) plot are due to effects of velocity crowding, i.e. to 
changing degrees of optical thinness, and a finite velocity dispersion of the 
gas. Furthermore, the line of sight kinematic point of symmetry with re- 
spect to the outer Galaxy is entirely coincident with the location of the sun. 
Thus, the most extreme velocities occur at galactic longitudes of I = 90° 
respectively I = 270° and show a drop-off in magnitude that is symmetrical 
with respect to the velocity minima at the center (I = 0°) and anticenter 
(I = 180°) directions. Note also, how in the 2 nd quadrant some emission is 
present at formally forbidden positive velocities, as is some emission in the 
3 rd at negative velocities, due to the velocity dispersion of nearby gas. 

5.2 Circular kinematics with warp 

Since the presence of a warp in the Galactic HI layer will alter the appearance 
of the velocity cutoff at b = 0° in those directions were the line of sight will 
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Figure 19: (l,v) plot at b = 0°, for the simulated data with circular kinematic 
symmetry, but a symmetrical warp with its maximum at 8 = 95° (upper 
panel), as well as a latitude integrated version for \b\ < 10° (lower panel), 
to show the effects this has on the influence of the warp. As in Figure 1, 
the peak intensity of the upper-panel display corresponds to a brightness 
temperature of 135 K and the lightest grey-scale corresponds to T;,=4 K. 
The lower-panel display was averaged over its b values and has a deviating 
grey-scale scaled between minimum and maximum intensities. 

leave the gas layer where it warps away maximally from the plane at radii 
smaller than those for which it reaches its extreme velocities, it is instructive 
to compare the observed (/, v ) plot at b = 0°, with an (/, v) plot of a model 
Galaxy described as in subsection 5.1.1, only now with a linear sinusoidal 
warp. The shape of the warp is such that a line through the midplane along 
a cylindrical wall through 8 and z at constant radius R is a perfect sine 
curve. The magnitude of the maximum of the warp starts at zero at R = 9 
kpc and increases from there with 110 pc for every kpc outwards, as was 
derived from the real data in section 4.2. For the simulations in this section 
we put the warp maximum at 9 = 95° . 

5.2.1 Conversion controlling 

In Figure 19 an (l,v) plot of the sort of Figure 18 was displayed, this time 
including a linearly rising sinusoidal warp. For comparison with a plot shown 
in Figure 23 in section 6.1 a latitude integrated version for \b\ < 10° was 
displayed as well. The upper panel shows indentations at the outer-Galaxy 
velocities in the heliocentric I directions where we expect the maximum warp 
along galactocentric 8. This complete Tb(l,b,v) data cube was processed in 
the same manner as the observed data cube to yield volume densities in 
cylindrical cuts through galactocentric coordinates of the sort displayed in 
the Figures 3-6. Since, in this case the input kinematics are, for certain, 
conveniently identical to those used in the velocity-to-distance processing, 
this should render a Galactic morphology that is consistent with the input 
model and, furthermore, draw our attention to any vagaries that can only 
have been inflicted by the method of converting the data cube and that can 
thus also be expected to blur our view of the Galaxy, displayed in the Figures 
3-6. The result of this conversion is displayed in galactocentric cylindrical 
cuts in Figure 20 for the constant radii R = 16 kpc and R = 26 kpc. 
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Figure 20: Mean HI volume densities in cylindrical cuts through the simu- 
lated data cube with circular kinematic symmetry, but a symmetrical warp 
with its maximum at 9 = 95° at galactocentric distances of R = 16 kpc (up- 
per) and 26 kpc (lower). The cylindrical-coordinate (R,9,z) data set was 
sampled at 1° intervals in 9 and 25-pc intervals in z; the velocity-to-distance 
transformation displayed here was circularly symmetric. The darkest grey- 
scale was scaled to the maximum densities at each of both radii, to be able 
to show details even at i?=26 kpc. Compare, for the exact derived densities, 
Figure 21 (middle panel). For reference a line tracing the galactic plane at 
6 = 0° was drawn. 



Figure 21: Vertical distance from the plane b = 0° of the Gauss-fitted peak 
HI density of the simulated data cube with circular kinematic symmetry, but 
a symmetrical warp with its maximum at 9 = 95° at a distance of R = 16 kpc 
(upper), to be compared with Figure 8, a 0-averaged (20° < 9 < 140°) HI 
midplane density with increasing radius (middle) (comparable with Figure 
12) and the flaring of the azimuthally averaged Gaussian width of the vertical 
HI density distribution (lower), comparable with Figure 15. The parameters 
were calculated in the same manner as their real-data counterparts, that is, 
after conversion of the data to galactocentric coordinates. 



Figure 22: Gaussian width (h) of the vertical HI density distribution (up- 
per) of the simulated data cube with circular kinematic symmetry, but a 
symmetrical warp with its maximum at 9 = 95° at a galactocentric distance 
of R = 12 kpc (comparable with Figure 14) and azimuthal distribution of 
HI midplane densities, at R = 10 kpc (lower) (comparable with Figure 10). 
Note that these parameters were calculated in the same manner as their 
real-data counterparts, that is, after conversion of the data to galactocen- 
tric coordinates. 
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5.2.2 Resulting vagaries 

A set of morphological parameters similar to those in section 4.1 were com- 
puted using the same procedures. They are displayed in Figures 21 and 
22. Several remarks can be made about the appearance of the cylindrical 
cuts in Figure 20 and on these parameters. First, it can be seen that the 
orientation of the warp is correct and also, at closer inspection, that the 
line of nodes is exactly at 9 = 5° and 9 = 185° (see also the upper panel of 
Figure 21). In fact, it was tested how the conversion analytics would deal 
with a warp maximum at exactly 9 = 0°, in view of the apparent coinci- 
dence of having the sun-center line so close to the observed line of nodes. 
It behaved properly and retrieved the exact input situation. Secondly, the 
gas layer has generally a constant thickness along azimuth and shows flaring 
from the R = 16 to the R = 26 plot, as required by the input. A non-input 
feature is clearly the density indentations at the center and anticenter di- 
rections, due to a mathematical singularity in the conversion analytics at 
9 = 0°, 180°. Also, the striking patches, stretched greatly in z, are caused 
by the processing, rather than inherent to the model. Its origin is of local 
gas that reaches velocities, due to its finite velocity dispersion, that corre- 
spond kinematically with radii in the far outer Galaxy. The densities are 
then derived from the received intensities according to the multiplefold dis- 
tances. It occurs only at the cardinal directions, as it is there where the 
line of sight velocities remain quite small, even at large radii. Note that 
in the case of kinematic symmetry, as is considered here, the patch is also 
symmetrically oriented around both directions, whereas in the Figures 3-6 
it was displaced. Compare, in this context, for instance, the upper panel of 
Figure 22 with that of Figure 14. That the patch reaches greater vertical 
heights above the plane in the center direction can be easily understood in 
terms of the large kinematic distance that local gas emission at any par- 
ticular latitude corresponds to, with its corresponding large z height, since 
it needs to traverse the sun-center radius twice first. Again, in the center 
direction, and in that of the anticenter to a lesser extent as well, a vagary 
is observed, at R = 16 kpc, displaying a progressive lower volume density 
when approaching 9 = 0° (and 9 = 180°). This effect is even more clear in 
the lower panel of Figure 22 and is attributable to a non-negligible optical 
depth due to velocity crowding. This concentrates all emission in a small 
range of frequencies and thus causes the spectra to saturate, resulting in an 
underestimation of the true volume density in spatial distribution. Finally, 
this under abundance near the cardinal directions changes into an overabun- 
dance for plots at very large radii; see the cylindrical cut at R = 26 kpc. 
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By the drop in density according to the input scale length, r s , at these dis- 
tances its magnitude remains only a fraction of the inner constant density 
disk value, as is illustrated in the middle panel plot of Figure 21. The con- 
tribution of emission far into the velocity dispersion wings from gas nearer 
to the Galactic center becomes non-trivial at such low densities. Generally, 
the effects induced by velocity crowding and velocity dispersion are stronger 
in the center direction than in anticenter direction, because there the inner 
solar radius gas emission confers at low positive and low negative velocities 
with the outer-Galaxy emission through its own finite velocity dispersion. 
From the middle panel of Figure 21, it can be seen that the input density of 
n Hi = 0.3 cm" 3 yields realistic volume density values after conversion. The 
lower panel in that figure was shown for the remarkably steep rising scale- 
height for smaller radii, something that is also observed in the analysis of 
the real data in this project, as well as in that of others. It is, obviously, 
not according to the original input, but due to a strong broadening of the 
elongated patches mentioned before when approaching Ro, thus biasing the 
azimuthally averaged scale height to greater values. 



6 Deviations from circular symmetry 

There is substantial long-established evidence for bilateral asymmetry in the 
outer-Galaxy gas layer. The kinematic cutoff at b = 0° corresponding to the 
far outskirts of the Galaxy occurs at velocities some 25 km /s more extreme 
in the southern-hemisphere data than in the northern (Burton, 1973 5 ). 
Additionally, the tightest constraint of radial velocities occurs 6° removed 
from the cardinal direction I = 180°; for I = 90°, circular symmetry would 
formally forbid any positive velocities (except those due to velocity disper- 
sion), while, in fact, the emission peaks at v = +6 km/s (Burton, 1972 4 ). 
Early on, Kerr (1962) suggested that the apparent difference between 
the inner-Galaxy northern and southern rotation curves derived from ter- 
minal velocity data could be due to an outward motion of the adopted local 
standard of rest of about 7 km/s. Blitz &; Spergel (1991) [2| reproduced the 
shape of the difference between velocity contours in the 0° < / < 180° and 
180° < I < 360° regions corresponding to a 1 K isotherm by adopting a 14 
km/s radial outward LSR motion and, assuming this to be kinematically 
equivalent to a more naturally explained radial outward component at the 
sun's location along global elliptical streamlines, they put forward a model 
of gas streamlines with decreasing ellipticity outwards, the outermost gas 
moving in almost circular orbits, with an angle a between the sun-center 
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line and the semi-major axis of the gas distribution of ~ 45° ± 20°, and an 
ellipticity e ~ 0.02, which they found to be consistent with some other ob- 
servables they examined (e.g. stellar kinematics, molecular gas). Kuijken & 
Tremaine (1991) and Kuijken (1992) 16 argued that from their exam- 
ination of the kinematics of the solar neighborhood they found little or no 
evidence for large-scale nonaxisymmetry. They note, however, that it would 
be very difficult to detect any ellipticity, if the sun happens to reside near 
a symmetry axis. These and other investigations were reviewed by Burton 
(1991) 0. 

The dynamical situation responsible for the observed asymmetries re- 
mains unclear, as does the extent to which the asymmetries might be spatial 
rather than kinematic. However, it is clear that even small deviations from 
kinematic symmetry are likely to have an observable impact on the HI pro- 
files, since their shape is very sensitive to small variations in the streaming 
motions (Burton, 1972) [3]. We will consider in what follows consequences 
of kinematical asymmetry only, although we note that it seems unlikely that 
this would be the exclusive cause of the observational situation in view of 
the common occurrence of true lopsidedness in other nearby spiral galaxies 
(Baldwin et al., 1980 1 j: Richter & Sancisi, 1994 |U). In any case, it is 
important to realize that many kinematic asymmetries are necessarily joined 
by a spatial counterpart (e.g. non-zero ellipticity, with the sun-center line 
anywhere in between the semi-major and semi-minor axes of the ellipsoid) 
and that a true structural lopsidedness, due to an m = 1 harmonic term in 
the Galactic potential, is always revealed by a perturbation in the line of 
sight velocity field (Swaters et al., 1999 EH]). 

We do not go into a full detailed treatment in order to pinpoint a definite 
Galactic morphology and its kinematics that can account for the observed 
deviations from circular symmetry, but instead we want to make use of our 
model simulations as a controlled input, showing the observational impli- 
cations of assumed real galactic asymmetries that have been suggested as 
explanations of the observed asymmetries on several aspects of observed 
spectra at the same time. In fact, this gives us a powerful instrument for 
putting constraints on various solutions to the asymmetry problem in the 
literature. 

In order to do so, we will first distinguish a set of observed asymmetries 
and establish a method to measure these asymmetries quantitatively. Then, 
after having found the real data's values of these parameters, we will adopt 
a list of possible morphologies and kinematics and measure the same set 
of parameters in the data as simulated according to these models. The 
measurements take place in an (/, v) -displayed set of spectra at b = 0°. 
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Figure 23: (l,v) plot at b = 0° (upper panel) as well as a latitude integrated 
version for |6| < 10° (lower panel), which can be considered free of any 
influence by deviations from the galactic plane. Like in Figure 1, the peak 
intensity of the upper panel corresponds to a brightness temperature around 
135 K and the lightest grey-scale corresponds to Tf,=4 K. The lower panel 
was averaged over its b values and has a deviating grey-scale scaled between 
minimum and maximum intensities. 

Some measurements shall be taken as well from a 6-integrated (I, v) display, 
for \b\ < 10°, to remove most of the warp's influence on these displays as is 
demonstrated in the next section. 

6.1 Asymmetries in the data 

We define spectral asymmetry as any feature or property of the T\, spectral 
(l,v) images, insofar as not corresponding to those according to simple pla- 
nar circularly symmetric rotation, with the assumed flat rotation curve as 
described in section 4. 

In the (I, b, v) data we choose to monitor the following set of parameters, 
when we reproduce galactic spectra according to different models. Most of 
these have been used before in investigations on Galactic asymmetry and 
some have been added by us, not always as an asymmetry, but also to put 
restrictions on the various explanations of observed asymmetry. We explain 
how we proceed to measure these phenomena and give the according values 
for the real data. 

1. The kinematic or spatial cutoff appears to occur at more extreme 
velocities in the part of the hemisphere accessible from the south than 
that from the north. 

2. The maximum velocities attained in the north and the south are not 
centered on I = 90°, 270° respectively, as would be expected in the 
symmetrical case. 

3. The shapes of the extreme velocity curves in the north and the south 
differ. 

4. The tightest velocity constraint towards the anticenter is displaced by 
a few km/s into the third quadrant. 
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5. The centroid or middle of this velocity distribution still appears cen- 
tered on I = 180°. 

6. The tightest velocity constraint towards the center is not displaced as 
it is towards the anticenter. 

7. The centroid of the velocity distribution around I = 0° is significantly 
displaced to positive velocities. 

8. Emission extends to higher positive velocities in the second quadrant 
than to corresponding negative velocities in the third. 

9. Finally, it appears that the velocity integrated column depth is large 
in the fourth quadrant with respect to the first, and that of the second 
quadrant is large compared to the third. 

We can further make remarks on the apparent asymmetries, as observed 
in the real data after conversion under the assumption of circular symmetry, 
in the (R, 6, z) frame, as can be seen in the Figures 3 through 6. They were 
not measured quantitatively. 

1. In general the South is observed to be endowed with a progressively 
higher abundance of gas with increasing radii than the North. 

2. An elongated patch of gas in the center direction, reaching high z, is 
observed to be displaced from I = 0° into the fourth quadrant. 

These effects can be real spatial anomalies (which would for point 2 be 
highly unlikely) or the result of the conversion itself under the, at least in 
part, wrong assumptions of circular symmetry. 

In particular the effect mentioned in point 1 could be reproduced if the 
kinematic situation would be such that gas on a certain radius from the 
Galactic center would give rise to a higher velocity along the line of sight in 
the south than in the north, be this due to its own velocity or due to that 
of the LSR 1 . 

Furthermore, the effect mentioned in point 2, when centered symmet- 
rically on I = 0°, is caused by local velocity-dispersed gas, entering at the 
center direction into negative velocities for the first quadrant and into posi- 
tive velocities for the fourth and hence accounted for originating from beyond 
the solar circle near I = 0° and scaled to densities accordingly (see also the 

1 In the Hartmann & Burton Leiden/Dwingeloo Survey the LSR was defined in terms 
of the Standard Solar Motion of 20 km/s toward (a, S) 1900 = (18\ 30°). 
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comments in the caption of Figure 13). This effect can also be seen in the 
simulation of section 5.2. The displacement can be easily conceived of by 
assuming a line of sight net positive velocity between the LSR and the local 
gas in the center direction, be this due to local streaming or to a radial 
velocity gradient by a non-circular global streaming. 

To quantify all these conditions we measure parameters as follows. For 
the items 4, 5, 6, and 7, we make use of the IDL routine GAUSSFIT as 
described in section 4.1 to fit a Gaussian to the velocity distribution in 
the (l,v) plot at certain longitudes around the cardinal directions I = 0° 
respectively I = 180°, after we have removed spectral features that are less 
than a factor e -2 times the maximum value at that specific longitude, 
that otherwise might have contaminated the outcome. The tightest velocity 
constraints are then taken to be the minima along a plot of the Gaussian 
widths varying with longitude, while the centroid of the velocity distribution 
simply is the location of the center of the Gaussian with longitude. For the 
items 1, 2, and 3, we only make use of the contour in the (l,v) display at 
Tj = 2 K. Item 1 is the maximum velocity along this contour over a broad 
longitude range around I = 270°, respectively the minimum around I = 90°. 
The longitudes where these extremes were recorded represent item 2. A 
parameter, Vdiff, is calculated from the difference between the magnitude 
of the negative velocities at I, for 0° < / < 180°, and the positive velocities 
at their co-longitudes, at 360° — /. This parameter corresponds to the shape 
of the extreme velocity curves, mentioned in item 3, and is positive for 
larger velocities in the south than in the north at any particular longitude. 
It is the same parameter that was used by Blitz & Spergel (1991) as 
the principal indication of their suggested outwardly decreasing ellipticity 
model. The emission mentioned in item 8 is probably due to local gas that, 
in a purely circularly symmetrical case, by its velocity dispersion would enter 
the (for its circular kinematics) forbidden quadrants in equal amounts (as 
was shown by the simulations of section five), but now has an extra positive 
line of sight velocity. If the emission were due to gas from larger radii in 
the second quadrant, shifted sufficiently to positive velocities to show up at 
v > km/s, it might surprise why this large-scale feature does not extend 
into the 1 st and the 3 rd quadrants, but instead disappears at their boundary 
quite abruptly. Also, it is shown below in section seven that a net radial 
inward (to the Galactic center) velocity of the LSR would reproduce the 
desired effect, but this is strongly ruled out by other (global) kinematic 
measures. Furthermore, a tangential error in the LSR would either increase 
the extent to which the emission enters the formally forbidden velocities 
for both quadrants or decrease them at the same time. Finally, item 9 is 
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Figure 24: Compilation of parameter plots derived from an (I, v) display at 
6 = 0°. The upper-left and middle plots show the net radial velocity of the 
peak 21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (<j) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity of 
a 2 K contour of longitudes versus their co-longitudes (360° — /), such that 
it is positive for larger velocities in the southern hemisphere. Finally, the 
lower-right plot is a measure of the percentage of excess column depth of the 
4 th quadrant with respect to the I s *, and of the 3 rd with respect to the 2 nd . 
In addition to this it followed that the most extreme velocity in the first two 
quadrants along a contour of 2 K in brightness temperature was —126 km/s 
occurring at / = 117° (27° displaced from / = 90°), while for the last two it 
was +144 km/s occurring at I = 281° (11° displaced from I = 270°), which 
renders their absolute ratio to be 1.1. Furthermore, emission from local gas 
enters along the same contour the formally forbidden second quadrant at an 
average velocity of 16.5 km/s, while in the third at —11.0 km/s, establishing 
an absolute ratio of 1.5. 

quantified by integrating over velocities with v < km/s for 0° < I < 180° 
and with v > km/s for their co-longitudes, at 360° — I. The / values 
are then subtracted from the (360° — /) values and this is divided by the 
(360° — I) values and then multiplied by 100 to give change in column depth 
in percentages. The integration cutoff at v = km/s does not significantly 
affect this parameter, N^iff, but instead the plot remains quite similar if we 
take the cutoff at v = ±10, ±20 or ±30 km/s, as was also shown by Blitz & 
Spergel. 

Figure 24 shows a compilation of parameter plots derived from an (l,v) 
display at b = 0° of the composite data cube and Figure 25 shows the same 
set of parameters for a latitude integrated (I, v) data set in the range \b\ < 
10°. Since the data grid points are separated by half a degree in latitude, this 
last measurement includes 41 cuts through the data of the sort that was used 
in the first. The change in the plots from b = 0° to |6| < 10° was monitored 
for the integrations \b\ < 1°, \b\ < 2°, and |6| < 5° as well and found to 
behave continuously and to settle fast towards the situation of the |6| < 10° 
plots. It can thus be relied upon that this plot reflects the global kinematic 
(and spatial) situation more accurately than any of the others, moreover so 
because effects due to the deviation of the HI layer from the plane b = 0° at 
larger radii, are being prevented, since emission at intensities large enough to 
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Figure 25: Same plots from an (I, u)-display at \b\ < 10° as in Figure 24 
are shown here. Again the upper-left and middle plots show the net radial 
velocity of the peak 21 cm emission with respect to the LSR towards the 
center respectively anticenter. The lower-left and middle plots show the 
width (a) of the total spectral distribution along any galactic longitude in 
the center respectively anticenter direction. The upper-right plot shows the 
difference in velocity of a 2 K contour of longitudes versus their co-longitudes 
(360° — I), such that it is positive for larger velocities in the southern hemi- 
sphere. Finally, the lower-right plot is a measure of the percentage of excess 
column depth of the 4 th quadrant with respect to the 1 st , and of the 3 rd with 
respect to the 2 nd . Corresponding values for the extreme velocity cutoffs at 
T b =2 K are -123 km/s at I = 105° (15° displaced from I = 90°) and 132 
km/s at / = 277° (7° displaced from / = 270°), with a ratio of 1.1. The local 
gas has averages of 15.2 km/s respectively —12.2, giving 1.2 as an absolute 
ratio. 

have an influence on the parameters under consideration here predominantly 
originates from within the covered latitude range. The difference between 
the plots of the centroid around I = 0° can therefore be attributed to a local 
irregularity in global streaming in the 6 = 0° case, something that is, in 
center direction, not implausible. It is reassuring for this reasoning that in 
the direction I = 180°, where great kinematic turmoil is not observed, the 
velocity centroid of both distributions have almost identical positions. The 
greater spectral width of the plot of tightest ^-constraint with longitude 
around / = 0° in the case of the 6 = 0° data cut with respect to the |6| < 10° 
version, while in the direction of longitudes around / = 180° they keep 
identical values, is a feature that appears to be regular and inherent to the 
dynamical situation and is reproduced in all of the simulations we made 
in this project. In particular, note the difference between Figure 27 and 
Figure 28, in this context. The change of the parameter vjiff along 0° < 
\l\ < 180° indicates that the outer-Galaxy line of sight radial velocities, 
measured at T^=2 K, are more extreme in the 4 th quadrant, and well into 
the 3 rd , with respect to those from the 1 st quadrant, and into the 2 nd . At 
certain directions (/ ~ 110° and I ~ 250°) this reverses and the velocities in 
the 2 nd quadrant become more extreme than those in the 3 rd . The column 
depth (Ndiff) varies similarly. Here, around I = 100°, there is a high relative 
peak to be seen of velocity integrated densities in the southern data with 
respect to the northern, a structure that has attributed to the Carina spiral 
arm by some authors, as was mentioned by Blitz & Spergel. In any case, this 
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Figure 26: Same plots from an (l,v) display at b = 0° as in Figure 24 are 
shown here for the simulated data according to a circularly symmetric model. 
Again the upper-left and middle plots show the net radial velocity of the 
peak 21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (a) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity of 
a 2 K contour of longitudes versus their co-longitudes (360° — I), such that 
it is positive for larger velocities in the southern hemisphere. Finally, the 
lower-right plot is a measure of the percentage of excess column depth of 
the 4 th quadrant with respect to the I s *, and of the 3 rd with respect to the 
2 nd . Corresponding values for the extreme velocity cutoffs at T^=2 K are 
-156 km/s at I = 90° (0° displaced from I = 90°) and 156 km/s at I = 270° 
(0° displaced from / = 270°), with a ratio of 1.0. The local gas has averages 
of 14.1 km/s respectively —14.1, giving 1.0 as an absolute ratio. 

sort of large structures seem to suggest that we should be cautious, when we 
try to explain the shapes of these parameters by global kinematic measures 
only, omitting other large scale influences like e.g. spiral structure. 

6.2 Symmetrical reference data 

As a reference to what would be the observational situation in case of per- 
fect circular symmetry or circular kinematic symmetry, but with a linear 
sinusoidal warp, and as yet another check on the behavior of the synthetic 
spectra, we performed the same asymmetry analysis as above on the con- 
trolled models that were described in section five. When considering perfect 
circular symmetry, there are indeed no deviations detected of the sorts cat- 
aloged in the previous subsection. This can be verified in Figure 26. This 
situation changes slightly though, with respect to the parameters v^iff and 
Ndiff, for the model including a linear sinusoidal warp (Figure 27). When 
we integrate the (l,v) display for \b\ < 10°, however, this difference almost 
totally vanishes, since then the warped layer is included as well (Figure 28). 
Note, as was mentioned above, how this also brings the two different tightest 
velocity constraints closer together. 
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Figure 27: Same plots from an (I, v) display at b = 0° as in Figure 24 are 
shown here for the simulated data according to a circularly kinematically 
symmetric model with a linear sinusoidal warp, with its maximum at 9 = 
95°. Again the upper-left and middle plots show the net radial velocity of the 
peak 21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (a) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity 
of a 2 K contour of longitudes versus their co-longitudes (360° — I), such 
that it is positive for larger velocities in the southern hemisphere. Finally, 
the lower-right plot is a measure of the percentage of excess column depth 
of the 4 th quadrant with respect to the I s *, and of the 3 rd with respect to 
the 2 nd . Corresponding values for the extreme velocity cutoffs at Tb=2 K 
are -133 km/s at / = 98° (8° displaced from / = 90°) and 136 km/s at 
I = 259° (-11° displaced from / = 270°), with a ratio of 1.0. The local 
gas has averages of 14.1 km/s respectively —14.1, giving 1.0 as an absolute 
ratio. 



Figure 28: Same plots from an (/, v) display at \b\ < 10° as in Figure 24 are 
shown here for the simulated data according to a circularly kinematically 
symmetric model with a linear sinusoidal warp, with its maximum at 9 = 
95° . Again the upper-left and middle plots show the net radial velocity of the 
peak 21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (<j) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity of 
a 2 K contour of longitudes versus their co-longitudes (360° — /), such that 
it is positive for larger velocities in the southern hemisphere. Finally, the 
lower-right plot is a measure of the percentage of excess column depth of the 
4 th quadrant with respect to the 1 st , and of the 3 rd with respect to the 2 nd . 
Corresponding values for the extreme velocity cutoffs at T^=2 K are —135 
km/s at I = 95° (5° displaced from / = 90°) and 135 km/s at I = 265° (-5° 
displaced from / = 270°), with a ratio of 1.0. The local gas has averages of 
14.1 km/s respectively —14.1, giving 1.0 as an absolute ratio. 
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6.3 Modeled asymmetries 

In general we can distinguish some categories of asymmetries with respect 
to deviations from simple planar circularly symmetric rotation, with the 
assumed flat rotation curve as described above. 

• Deviations from planar rotation. 

• Deviations from a flat rotation curve. 

• Deviations in the LSR. 

• Deviations from circular streaming. 

The only deviation with immediate observational consequences that is 
quite firmly established is the one that is first mentioned and has for a 
symmetrical warp been tested already in the previous subsection. The pos- 
sibilities that were further investigated here were the following: 

• Asymmetrical warp 

• Radial velocity LSR inward/outward from the GC. 

• Global circular streaming, with sun on elliptical movement. 

• Global elliptical streaming. 

• Small radial velocity of local gas in the solar neighborhood. 

• Linear change in ellipticity from solar radius outward. 

With respect to deviations from a flat rotation curve, we can comment 
as follows. A linear or monotonically and smoothly varying rotation curve 
will cause only scaling effects due to the different velocity-to-distance trans- 
formation. Local perturbations with R, however, or, for that matter, with 
6, certainly will cause the appearance of the spectra and, after conversion, 
that of the derived Galactic morphology, to change significantly. Since we 
consider global parameters only here, we did not perform such simulations. 
During our investigations various different values for the parameters that 
occur in the models were tried out to monitor the response of the generated 
spectra to these changes, as will be described in the next section. 
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7 Observational implications of model asymme- 
tries 

After having cataloged, in the previous section, a set of observed asymme- 
tries and constraints and listed a number of morphologies and kinematics 
that could mimic these asymmetries, we will proceed to pursue the obser- 
vational implications of certain specific models for the observed spectra. It 
has been shown which properties of the spectra require that the assumptions 
of circular spatial and kinematic symmetry be abandoned. In this section 
it will become apparent which modeled asymmetries produce spectra that 
point in the right direction. For clearness, we will monitor one change with 
respect to the standard model of section 5.1 at a time, except in 7.7, when 
we consider global elliptical streaming in combination with an asymmetrical 
warp. This, to give galactocentric cylindrical plots that are comparable with 
the real data, since the asymmetrical warp is quite firmly observationally 
established (and it is the only one in this section to have been so). We are 
putting limits on a possible error in the determination of the motion of the 
Local Standard of Rest as well as on possible ellipticity in combination with 
different orientations of the semi-major axis with respect to the line con- 
necting the sun with the Galactic center. Such galactic asymmetries could 
mimic some of the observed bilateral asymmetries in the (l,b,v) data. We 
conclude this section with a consideration of the effects of a velocity-to- 
distance conversion, using circular kinematics, on the (I, b, v) spectra of a 
Galaxy that, in fact, rotates along global elliptical streamlines. 

7.1 Circular symmetry with asymmetrical warp 

The first deviation from the model in section 5.1 is that of an asymmetrical 
warp. It should be noted that this deviation is less hypothetical than the 
others that will follow, since it has been accurately measured. Even though 
this measurement was performed in the data cube that was constructed 
using circular kinematics instead of some possibly better, but yet unknown, 
alternative, it cannot be expected to differ more than of the order of, at 
most, ten percent from reality. The shape of the warp is again sinusoidal, 
but this time the amplitude differs in the part of the hemisphere that is 
accessible from the north from that of the south, according to the curves 
that are plotted in Figure 16. The maximum of displacement from the 
plane 6 = 0° was set to be attained at = 120°. It was separated from the 
minimum by 180°. 

An (l,v) plot of the resulting spectra is shown in Figure 29. It is clearly 
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Figure 29: Composite spectral image at b = 0°, for a model with circular 
kinematical symmetry, but asymmetrical warp with its maximum at 8 = 
120°, as described in the text. As in Figure 1, the peak intensity of the 
upper panel display corresponds to a brightness temperature of 135 K and 
the lightest grey-scale corresponds to T^=4 K. 

seen that in the 3 rd and 4 th quadrants, by the folding-back property of the 
asymmetrical warp in the south, velocities are recorded along the lowest 
brightness temperature contours that are more extreme than those of the 
1 st and 2 nd quadrants. Furthermore, in the compilation of parameter plots, 
that check for asymmetry, in Figure 30, it is interesting to note the large 
impact this has on the parameters vjiff and N^ff, and on the ratio of the 
extreme velocity cutoffs, for a cross-cut at b = 0°. Taking into account 
the relative certainty with which we should accept the asymmetrical warp 
effects, they should be imagined to reside implicitly in the same real data's 
6 = 0° parameter plots. However, it was demonstrated earlier (in section 
6.2) that for a |6| < 10° plot these effects largely disappear. Also, mind that 
the distortions of the above mentioned parameters from symmetry due to 
the asymmetrical warp will have to be superimposed on those resulting from 
the following investigations, if comparison is to be made with the plots in 
Figure 24. Particularly, it is encouraging that the N^iff parameter in this 
figure is quite indented for the first few tens of degrees in |/|, when compared 
with its |6| < 10° counterpart, consistent with the curve below 0% of N^iff 
in Figure 30. To visualize its asymmetrical character, the warp was also 
displayed in a cylindrical cut through z and 9 at the galactocentric distance 
of R = 20 kpc. 

7.2 Radial velocity LSR 

Superimposed upon the standard flat plane-parallel simulated Galaxy with 
circular kinematic and spatial symmetry (section 5.1) is the effect of a 10 
km/s motion of the LSR radially outward from the Galactic center. Figure 
32 shows the spectral situation of this model, again for 6 = 0°. From this 
figure and from the compilation of parameter plots in Figure 33 it can clearly 
be seen that both in the center and in the anticenter directions the centroid 
of the spectrum has an off-set exactly as great as the initial radial velocity in 
absolute magnitude. The Vdiff parameter shows a reasonable fit to the real 
data (Figure 25), although a slightly better fit would be attained for a higher 
radial outward velocity. This is exactly what Blitz & Spergel derived for the 
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Figure 30: Same plots for b = 0° as in Figure 24, for the case of circular 
kinematical symmetry with an asymmetrical warp with a maximum at 9 = 
120°. Again the upper-left and middle plots show the net radial velocity 
of the peak 21 cm emission with respect to the LSR towards the center 
respectively anticenter. The lower-left and middle plots show the width (<r) 
of the total spectral distribution along any galactic longitude in the center 
respectively anticenter direction. The upper-right plot shows the difference 
in velocity of a 2 K contour of longitudes versus their co-longitudes (360° — 
I), such that it is positive for larger velocities in the southern hemisphere. 
Finally, the lower-right plot is a measure of the percentage of excess column 
depth of the 4 th quadrant with respect to the 1 st , and of the 3 rd with respect 
to the 2 nd . Corresponding values for the extreme velocity cutoffs at T},=2 
K are -111 km/s at / = 89° (-1° displaced from I = 90°) and 156 km/s at 
/ = 270° (0° displaced from I = 270°), with a ratio of 1.4. The local gas has 
averages of 14.1 km/s respectively —14.1, giving 1.0 as an absolute ratio. 



Figure 31: Mean HI volume densities in a cylindrical cut through the simu- 
lated data cube at a galactocentric distance of R = 20 kpc, for a model with 
circular kinematical symmetry, but asymmetrical warp as described in the 
text. The warp maximum here was set at = 120°. From the intersection 
of the Galactic midplane with the plane established by b = 0° it can be seen 
that the line of nodes resulting from this warp maximum in the case of a 
sinusoidal warp is not consistent with what is seen in the real data (Figures 
2-5). 



7 OBSERVATIONAL IMPLICATIONS OF MODEL ASYMMETRIES^ 



Figure 32: Composite spectral image at b = 0°, for a model with circular 
symmetry but a radial outward velocity of 10 km/s as described in the text. 
As in Figure 1, the peak intensity corresponds to a brightness temperature 
of 135 K and the lightest grey-scale corresponds to T{,=4 K. 

Figure 33: Same plots for b = 0° as in Figure 24, for the case of circular 
symmetry with a radial outward velocity of 10 km/s. Again the upper-left 
and middle plots show the net radial velocity of the peak 21 cm emission with 
respect to the LSR towards the center respectively anticenter. The lower- 
left and middle plots show the width (a) of the total spectral distribution 
along any galactic longitude in the center respectively anticenter direction. 
The upper-right plot shows the difference in velocity of a 2 K contour of 
longitudes versus their co-longitudes (360° — /), such that it is positive for 
larger velocities in the southern hemisphere. Finally, the lower-right plot is 
a measure of the percentage of excess column depth of the 4 th quadrant with 
respect to the 1 st , and of the 3 rd with respect to the 2 nd . Corresponding 
values for the extreme velocity cutoffs at T{,=2 K are —157 km/s at I = 93° 
(3° displaced from I = 90°) and 157 km/s at I = 273° (3° displaced from 
/ = 270°), with a ratio of 1.0. The local gas has averages of 7.4 km/s 
respectively —20.9, giving 0.4 as an absolute ratio. 

same parameter in their simulations. As a best fit they found v ra d= 14 km/s 
outward. However, our other parameters strongly oppose such an outcome. 
For, not only would this increase the velocity of the centroid with respect 
to the center direction, which is not desirable in view of Figure 25, but also, 
and this is most important, does it contradict the observational evidence of 
little or no velocity with respect to the anticenter of the spectral centroid 
found there. On the other hand it is possible to reconcile N^ff with the 
observations, but it will be shown below that this is also possible for other 
models. The same is true for v^iff- Finally, by itself not conclusive, but 
still a matter of concern, is that the ratio of the average local gas velocities 
shows a directly opposite movement to what is observed. Note, that in case 
of an inward velocity all effects would be the same, but reversed. Through 
the above mentioned reasoning we rule out any correction to the LSR used 
in the Leiden/Dwingeloo survey greater than a few km/s. 
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Figure 34: Construction of the rotational situation used in the solar elliptical 
movement simulation. For clarity, the ellipticity is greatly exaggerated. 

Figure 35: Composite spectral image at b = 0°, for a model with circular 
kinematics, but with the sun moving on a slightly elliptical orbit with e = 
0.05 and a = 45°, as described in the text. As in Figure 1, the peak intensity 
of the upper panel display corresponds to a brightness temperature of 135 
K and the lightest grey-scale corresponds to T^=4 K. 

7.3 Solar elliptical movement 

In the case of a solar elliptical movement, but keeping the remainder of the 
Galaxy identical to that in section 5.1, the effects mimic in part those for 
a solar radial velocity. The rotational situation is shown in Figure 34 for a 
face-on view of the Galaxy as seen from the north pole. In Figures 35 and 
36 the results are shown for the case of an ellipticity of e=0.05 and an angle 
between the sun-center line and the semi-major axis of its orbit of a = 45°. 
The line of sight component of the sun's orbital velocity immediately creates 
a non-acceptable shift with respect to the anticenter centroid and even a 
double shift in the center direction. At the same time the ratio of the local 
gas average velocities again moves in the wrong direction. Any measure 
to reduce both shifts in order to get within reach of the real observables, 
reduces the effects on and N^ff correspondingly and thus seems not 
very fruitful. These measures include a decrease in ellipticity and a change of 
a towards both 0° and 90°. For these reasons, we rule out the solar elliptical 
movement according to the model as discussed as a natural explanation of 
the observed asymmetries. 

We must note that all alternatives that imply a motion of the LSR with 
respect to its immediate neighborhood are likely to be refuted by the same 
constraints mentioned above, since the bigger part of all the 21 cm intensity 
is emitted by gas that resides within a few kpc distance from the sun. The 
parameter v^i // on the other hand originates from the outer Galaxy only. As 
we shall see, is in the case of global elliptical streaming the velocity difference 
with respect to the immediate solar neighborhood limited, whereas with 
respect to the region outside the sun-center radius in the direction of I = 0° 
it is considerable. This is in direct correspondence with the observed data. 
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Figure 36: Same plots for b = 0° as in Figure 24, for the case of circular kine- 
matics with the sun moving on a slightly elliptical orbit, with e = 0.05 and 
a = 45°. Again the upper-left and middle plots show the net radial velocity 
of the peak 21 cm emission with respect to the LSR towards the center re- 
spectively anticenter. The lower-left and middle plots show the width (a) of 
the total spectral distribution along any galactic longitude in the center re- 
spectively anticenter direction. The upper-right plot shows the difference in 
velocity of a 2 K contour of longitudes versus their co-longitudes (360° — I), 
such that it is positive for larger velocities in the southern hemisphere. Fi- 
nally, the lower-right plot is a measure of the percentage of excess column 
depth of the 4 th quadrant with respect to the 1 st , and of the 3 rd with respect 
to the 2 nd . Corresponding values for the extreme velocity cutoffs at Tf,=2 
K are -161 km/s at I = 92° (2° displaced from I = 90°) and 151 km/s at 
/ = 272° (2° displaced from I = 270°), with a ratio of 0.9. The local gas has 
averages of 10.5 km/s respectively —17.6, giving 0.6 as an absolute ratio. 

Figure 37: Construction of the rotational situation used in the global ellipti- 
cal streaming simulation. For clearness, the ellipticity is greatly exaggerated. 

7.4 Global elliptical streaming 

The model incorporating global elliptical streaming is further clarified in 
Figure 37. Reviewed here is the situation for e=0.05 and a = 45°. Figures 
38 and 39 show the results of its analysis. It is immediately seen that 
this model combines positive results for all parameters. Thus, it causes 
the required positive velocity displacement for the centroid in the center 
direction, and still keeps that in the anticenter direction fixed at zero. Both 
Vdiff and Ndiff show a distortion that can, in combination with that of the 
asymmetrical warp, be united with the situation for the real data. On top, 
the ratio of the local gas average velocities agrees more here to that observed. 
We cannot give a fixed set of variables here that represent a best fit, but it 
is clear that this is the only one of our models that accounts for the general 
trends without contradicting others. In fact, it has not been found possible 
by any other model than global elliptical streaming, through the simulations 
we used, to reproduce the required positive velocity of the center direction 
centroid and at the same time keep the centroid in the anticenter direction 
fixed at v=0 km/s within a few km/s. Finally, if the angle a were to be 
set at negative values all plots would show the same curves, but inverted. 
It is interesting to see that the conclusion resulting from this subsection 
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Figure 38: Composite spectral image at b = 0°, for a model with global 
elliptical streaming as described in the text. As in Figure 1, the peak inten- 
sity of the upper panel display corresponds to a brightness temperature of 
135 K and the lightest grey-scale corresponds to T&=4 K. 

Figure 39: Same plots for b = 0° as in Figure 24, for the case of global 
elliptical streaming. Again the upper-left and middle plots show the net 
radial velocity of the peak 21 cm emission with respect to the LSR towards 
the center respectively anticenter. The lower-left and middle plots show the 
width (a) of the total spectral distribution along any galactic longitude in 
the center respectively anticenter direction. The upper-right plot shows the 
difference in velocity of a 2 K contour of longitudes versus their co-longitudes 
(360° — I), such that it is positive for larger velocities in the southern hemi- 
sphere. Finally, the lower-right plot is a measure of the percentage of excess 
column depth of the 4 th quadrant with respect to the I st , and of the 3 rd with 
respect to the 2 nd . Corresponding values for the extreme velocity cutoffs at 
T b =2 K are -166 km/s at / = 90° (0° displaced from / = 90°) and 147 km/s 
at / = 271° (1° displaced from I = 270°), with a ratio of 0.9. The local 
gas has averages of 14.1 km/s respectively —14.1, giving 1.0 as an absolute 
ratio. 

invokes a ellipticity and an orientation of the inflicted morphology, that is 
consistent with the model that Blitz & Spergel put forward. Although it 
can be seen that their assumption, that the kinematic influence on Vdiff of a 
radial component of the solar motion along elliptical streamlines is identical 
to a true radial velocity of the LSR, is not generally valid. 

7.5 Local streaming 

Since there was no global model found that could naturally explain the 
greater than one value of the absolute ratio of local gas average velocities 
without going against other strong constraints, a non-global explanation 
was thought of that could. This was that of local streaming. The model 
is again the Galaxy of section 5.1, only this time with an area around the 
sun's location within a radius of 250 pc that has a velocity radially outward 
from the Galactic center of 7 km/s. It is important to realize that by our 
parameters in the compilation of Figure 41 only the response of the data in 
the anticenter direction is monitored, and that, in fact, in the center direction 
the reasoning around the asymmetrical, elongated patch of gas (section 6.1) 
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Figure 40: Composite spectral image at b = 0°, for a model with local 
streaming as described in the text. As in Figure 1, the peak intensity cor- 
responds to a brightness temperature of 135 K and the lightest grey-scale 
corresponds to Tb=4 K. 

Figure 41: Same plots for b = 0° as in Figure 24, for the case of local stream- 
ing. Again the upper-left and middle plots show the net radial velocity of 
the peak 21 cm emission with respect to the LSR towards the center respec- 
tively anticenter. The lower-left and middle plots show the width (a) of 
the total spectral distribution along any galactic longitude in the center re- 
spectively anticenter direction. The upper-right plot shows the difference in 
velocity of a 2 K contour of longitudes versus their co-longitudes (360° — I), 
such that it is positive for larger velocities in the southern hemisphere. Fi- 
nally, the lower-right plot is a measure of the percentage of excess column 
depth of the 4 th quadrant with respect to the 1 st , and of the 3 rd with respect 
to the 2 nd . Corresponding values for the extreme velocity cutoffs at T^=2 
K are -156 km/s at / = 90° (0° displaced from / = 90°) and 156 km/s at 
/ = 270° (0° displaced from I = 270°), with a ratio of 1.0. The local gas has 
averages of 18.2 km/s respectively —11.0, giving 1.6 as an absolute ratio. 

would argue for a motion away from the sun. In Figures 40 and 41 it can be 
seen that local streaming reaches the right objectives and yet leaves the other 
observables almost unaffected. When we increase the radius of the local 
area that is given a net velocity, the parameters increasingly respond to it in 
accordance with an error in the LSR, as more overall HI intensity is taken to 
have a distorted velocity. As a local effect one would need to incorporate all 
local effects that are associated with the solar neighborhood on such scales, 
in order to derive valid conclusions about Galactic properties. As such it 
is not meant to describe the solar neighborhood in this section. However 
this model is taken to be suggestive of a remarkability in the analysis of the 
parameters we monitor, in view of an observed asymmetry. 

7.6 Decreasing versus constant ellipticity 

Finally, we comment upon the difference of constant ellipticity versus de- 
creasing ellipticity, as advocated by Blitz & Spergel. The same global ellip- 
tical model as in the subsection 7.5 above was altered to let the ellipticity 
decrease linearly beyond R =9 kpc to zero at R =25 kpc. A parameter 
compilation plot was made, as shown in Figure 42, comparable to that in 
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Figure 42: Same plots for b = 0° as in Figure 24, for the case of an outward 
decreasing elliptical streaming. Again the upper-left and middle plots show 
the net radial velocity of the peak 21 cm emission with respect to the LSR 
towards the center respectively anticenter. The lower-left and middle plots 
show the width (a) of the total spectral distribution along any galactic 
longitude in the center respectively anticenter direction. The upper-right 
plot shows the difference in velocity of a 2 K contour of longitudes versus 
their co-longitudes (360° — such that it is positive for larger velocities in 
the southern hemisphere. Finally, the lower-right plot is a measure of the 
percentage of excess column depth of the 4 th quadrant with respect to the 
I s *, and of the 3 rd with respect to the 2 nd . Corresponding values for the 
extreme velocity cutoffs at T^=2 K are —160 km/s at I = 92° (2° displaced 
from I = 90°) and 160 km/s at / = 272° (2° displaced from I = 270°), with 
a ratio of 1.00. The local gas has averages of 14.0 km/s respectively —14.2, 
giving 1.0 as an absolute ratio. 

Figure 39. The slight differences might be pleading in favor of the decreasing 
cllipticity, though not very strongly on the basis of our analysis. For com- 
parison, we tried a version with increasing ellipticity from the solar radius 
onwards, with negative results. 

7.7 Conversion remarkabilities 

There is a compelling argument that could make it worthwhile to test a 
conversion to galactocentric (R,9,z) coordinates from heliocentric (l,b,v) 
coordinates taking as an input model a combination of the best options of the 
alternatives that we went through in the subsections above. This would be 
a combination of the asymmetrical warp and the global elliptical streaming 
with a = 45° and an ellipticity of e=0.05. This argument is, of course, that if 
our Galaxy were to be of such morphology, with corresponding kinematics, 
then at the velocity-to-distance conversion, assuming the wrong type of 
kinematics, errors would be made with respect to the distances that are 
then attributed to certain emissivities and, what is at least as bad, volume 
densities are calculated accordingly. As in section 5.2, we thus constructed 
the complete (l,b,v) data cube for the combined model mentioned here. 
Resulting plots for spectra through v and / for 6 = 0° and for \b\ < 10° 
are shown in Figure 43. Compare the two compilation parameter plots in 
Figures 44 and 45 that were derived from them. It is instructive to note 
that the velocity of the centroid of the spectrum in the center direction is 
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Figure 43: (l,v) plot at b = 0°, for the simulated data with elliptical stream- 
ing and asymmetrical warp with its maximum at 6 = 95° (upper panel), 
as well as a latitude integrated version \b\ < 10° (lower panel), to show the 
effects this has on the influence of the warp. As in Figure 1, the peak inten- 
sity of the upper panel display corresponds to a brightness temperature of 
135 K and the lightest grey-scale corresponds to Tf,=4 K. The lower panel 
display was averaged over its b values and has a deviating grey-scale scaled 
between minimum and maximum intensities. 

lessened significantly when averaged over several (41 to be exact) cross-cuts 
in latitude. Compared to Figure 25, this ellipticity suddenly does not seem 
very large anymore, but instead quite appropriate. Also, note, how the 
effect of the asymmetrical warp is greatly reduced for the case \b\ < 10°. 
Now let us consider the remarkabilities concerning the velocity-to-distance 
data processing. As one can see in Figure 46, the discrepancy between the 
input model and its display in (R, 0, z) after conversion is by far greater 
than that, as described in section 5.2.2, of the case of identical model and 
conversion kinematics. Since, in addition to some of the vagaries that were 
commented upon there, we see some striking new ones. First, it appears 
that, just as in the Figure 3-6 displays, the 4 th quadrant is overabundant 
in gas with respect to the I s * and similarly the 2 nd quadrant displays a 
higher volume density than the 3 rd , especially for larger radii. Secondly, the 
rise to high z from the 2 nd quadrant to the 1 st of the HI gas layer and its 
subsequent detachment of that in the 4 th quadrant. A feature that, again 
more towards higher radii, is also seen in the real data. Finally, the apparent 
larger overall volume densities of the southern hemisphere data with respect 
to the northern. This is perceived to be due to the kinematic effect that was 
explained in section 6.1. 

8 Summary 

In this project the Leiden/Dwingeloo data was used as a main constituent in 
a composite data cube to display the outer reaches of our Galaxy. To account 
for any possible distortion of the warp description due to anomalous velocity 
gas the principle HVC complexes as identified by Wakker & van Woerden 
were masked out before the velocity-to-distance conversion took place. For 
this conversion we used the IAU standards of Ro=8.5 kpc and ©o=220 km/s. 
Due to the properties of the Leiden/Dwingeloo survey, at the same time a 
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Figure 44: Same plots from an (I, v) display at b = 0° as in Figure 24 
are shown here for the simulated data according to a model with elliptical 
streaming and asymmetrical warp with its maximum at 9 = 95°. Again 
the upper-left and middle plots show the net radial velocity of the peak 
21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (a) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity 
of a 2 K contour of longitudes versus their co-longitudes (360° — I), such 
that it is positive for larger velocities in the southern hemisphere. Finally, 
the lower-right plot is a measure of the percentage of excess column depth 
of the 4 th quadrant with respect to the I s *, and of the 3 rd with respect to 
the 2 nd . Corresponding values for the extreme velocity cutoffs at Tb=2 K 
are -133 km/s at / = 98° (8° displaced from / = 90°) and 136 km/s at 
I = 259° (-11° displaced from / = 270°), with a ratio of 1.0. The local 
gas has averages of 14.1 km/s respectively —14.1, giving 1.0 as an absolute 
ratio. 



Figure 45: Same plots from an (l,v) display at \b\ < 10° as in Figure 24 
are shown here for the simulated data according to a model with elliptical 
streaming and asymmetrical warp with its maximum at 9 = 95°. Again 
the upper-left and middle plots show the net radial velocity of the peak 
21 cm emission with respect to the LSR towards the center respectively 
anticenter. The lower-left and middle plots show the width (a) of the total 
spectral distribution along any galactic longitude in the center respectively 
anticenter direction. The upper-right plot shows the difference in velocity of 
a 2 K contour of longitudes versus their co-longitudes (360° — /), such that 
it is positive for larger velocities in the southern hemisphere. Finally, the 
lower-right plot is a measure of the percentage of excess column depth of the 
4 th quadrant with respect to the I s *, and of the 3 rd with respect to the 2 nd . 
Corresponding values for the extreme velocity cutoffs at T&=2 K are —135 
km/s at I = 95° (5° displaced from / = 90°) and 135 km/s at I = 265° (-5° 
displaced from / = 270°), with a ratio of 1.0. The local gas has averages of 
14.1 km/s respectively —14.1, giving 1.0 as an absolute ratio. 
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Figure 46: Mean HI volume densities in cylindrical cuts through the sim- 
ulated data cube with elliptical streaming and asymmetrical warp with its 
maximum at 9 = 95° at galactocentric distances of R = 16 kpc (upper) and 
26 kpc (lower). The cylindrical-coordinate (R, 6, z) data set was sampled at 
1° intervals in 9 and 25-pc intervals in z; the velocity-to-distance transfor- 
mation displayed here was circularly symmetric. The darkest grey-scale was 
scaled to the maximum densities at each of both radii, to be able to show 
details even at i?=26 kpc. For reference a line tracing the galactic plane at 
6 = 0° was drawn. The warp maximum here was set at 6 = 95°. 

more sensitive and more detailed description of the warp and its flaring na- 
ture was given than had been possible before. A range of warp parameters 
was further quantified and compared to previous investigations. In order 
to follow the method of converting from heliocentric (l,b,v) coordinates to 
galactocentric (R, 0, z) coordinates a complete set of synthetic spectra was 
constructed on the same grid as the Leiden/Dwingeloo survey and subse- 
quently processed in identical manner as the real data. Several observational 
vagaries that resulted from this test conversion have been identified in the 
real spectra. In addition controlled conversion has been done for a range 
of hypothetical Galaxy models to search for resulting remarkabilities that 
can be associated with the observational situation, since it has been long 
established that the spectra show asymmetries with respect to global circu- 
lar spatial and kinematic symmetry. In order to monitor this quantitatively, 
we choose a fixed set of parameters in the (l,b,v) data and plotted their 
measured values in a single compilation display. This was done for the real 
data as a reference as well. It appeared that many different models create 
an effective significant motion of the LSR, which is by our parameter plots 
strongly confined to very low velocities. Judging on a global scale these mod- 
els can therefore be ruled out. To the observational situation, both directly 
in heliocentric coordinates and in a more subtle (kinematical) way in galac- 
tocentric coordinates (after conversion using circular symmetry), remarks 
can be made that agree with those that follow from a modeled elliptical 
input Galaxy. This outer Galaxy ellipticity is consistent with that found by 
Blitz & Spergel. We note, furthermore, that it has the same orientation as 
the gas streamlines that resulted from the modeling of the motions in the 
inner few kpc by Liszt & Burton (1980) [IS] 
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